Cosmological Evolution of Interacting Phantom Energy with Dark Matter 
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We investigate the cosmological evolution of an interacting phantom energy model in which the 
phantom field has interaction with the dark matter. We discuss the existence and stability of scaling 
solutions for two types of specific interactions. One is motivated by the conformal transformation 
in string theory and the other is motivated by analogy with dissipation. In the former case, there 
exist no scaling solutions. In the latter case, there exist stable scaling solutions, which may give a 
phenomenological solution of the coincidence problem. Furthermore, the universe either accelerates 
forever or ends with a singularity, which is determined by not only the model parameters but also 


the initial velocity of the phantom field. 


PACS numbers: 98.80.Cq, 98.80.-k 


I. INTRODUCTION 


Scalar field plays an important role in modern cosmol- 
ogy. The dark energy can be attributed to the dynamics 
of a scalar field, for instance quintessence al , which 
convincingly realize the present accelerated expansion by 
using late-time attractor solutions, in which the scalar 
field mimics the perfect fluid in a wide range of parame- 
ters. Much attention has been drawn to the case of expo- 
nential potentials. The exponential potentials allow the 
possible existence of scaling solutions in which the scalar 
field energy density tracks that of the perfect fluid (so 
that at late times neither component can be negligible). 
In particular, a phase-plane analysis of the spatially flat 
FRW models showed that these solutions are the unique 
late-time attractors whenever they exist More- 
over, exponential potentials appear naturally in the low 
dimensional models of string/M-theory [a]. 


The recent SNe data seem to favor the _dark energy 
with the present equation of state w < —1 [dj To obtain 
w < —1, phantom field with a negative kinetic term may 
be a simplest implementing and can be regarded as one of 
interesting possibilities describing dark energy m. The 
physical background for phantom type of matter with 
strongly negative pressure would be found in string the- 
ory |8]. Phantom field may also arise in higher-order 
theories of gravity [9], Brans-Dicke and non-minimally 
coupled scalar field theories [o]. The cosmological mod- 
els which allow for phantom matter appear naturally in 
k-essence models mi In spite of the fact that the field 
theory of phantom fields encounters the problem of sta- 
bility which one could_try to bypass by assuming them 
to be effective fields [12], it is nevertheless interesting 
to study their cosmological implication. Recently, there 
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have been many relevant studies on this topic Ep 


The physical properties of phantom energy are rather 
weird, as they include violation of the dominant energy 
condition and increasing energy density with the expan- 
sion of the universe. The latter ultimately leads to an un- 
wanted future singularity called big rip. This singularity 
is characterized by the divergence of the scale factor in 
a finite time in future [14]. To avoid the cosmic dooms- 
day, some phantom field models were proposed E} It 
requires a special class of phantom field potentials with 
a local maximum. Moreover, the energy density of the 
phantom field increases with time, while the energy den- 
sity of the matter fluid decreases as the universe expands. 
Why are the energy density of dark matter and the phan- 
tom energy density of the same order just at the present 
epoch? This coincidence problem becomes more diffi- 
cult to solve in the phantom model. Therefore, from this 
point of view the cosmological scaling solution would be 
desirable for the history of the universe. Throughout this 
paper we use “scaling solution” as a meaning that the en- 
ergy densities of the phantom field and the dark matter 
are proportional. However, as shown in Ref. Ed, there 
exist no scaling solutions because the phantom energy 
increases while the matter energy decreases with time. 
But in the presence of the suitable interaction this case 
could be realized easily in Ref. Li. In this paper we in- 
vestigate the stability and existence of scaling solutions 
in the scenario of interacting phantom energy with dark 
matter. We consider two phenomenological models. One 
is motivated by the conformal transformation from the 
Jordan frame to Einstein frame in string theory and the 
other is motivated by analogy with dissipation. In the 
former model there exist no scaling solutions. However, 
in the latter model a phase-plane analysis shows that 
there exist two kinds of stable scaling solutions, which 
lead to two different fates of the universe. The universe 
either accelerates forever or ends with a big rip, which is 
determined by not only the model parameters but also 


the initial velocity of the phantom field. We will also 
discuss the physical consequences of these results. 


II. INTERACTING PHANTOM ENERGY WITH 
DARK MATTER 


Let us consider a universe model where both the phan- 
tom field ¢ and the dark matter pm are present. The 
Friedmann equation in a spatially flat FRW metric can 
be written as 

2 


H? = (pp + pm), (1) 


where K? = 87Gy is the gravitational coupling and the 
energy density and pressure, pp and P,, of the homoge- 
neous phantom field ¢ are given by 


Pp = -2+ VO) (2) 
P, = -5#-Vi9), (3) 


respectively, in which V(¢) is the phantom field poten- 
tial. We postulate that the two components, pp and pm, 
interact through the interaction term Q according to 


Pr + 3H (Pm + Pm) = Q, (4) 
bp + 3H (pp + Py) = -Q. (5) 


Suppose the dark matter possesses the equation of 
state P,, = 0. The dynamics of the phantom field with 
an exponential potential 


V ($) = Vo exp(—Ax¢) (6) 


has been analysed in Ref. ig. We assume the dimension- 
less constant A is positive since we can make it positive 
through ¢ — —@ if A < 0. Here we generalize the analy- 
sis to the case in which the phantom field has interaction 
with the dark matter. As we will see, this gives rise to 
some interesting novel features. Following Ref. |3], we 
define the following dimensionless variables 
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Notice that x”, y? and z? give the fraction of total en- 
ergy density carried by the field kinetic energy, the field 
potential energy and the dark matter, respectively. Thus 
the fractional densities of pp and pm can be written as 
Qp = —x? +y? and Qm = z?, respectively. The evolution 
equations (4) and (5) can be written as the following set 
of equations: 


(7) 
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where the prime denotes a derivative with respect to the 
logarithm of the scale factor, N = Ina, and the Fride- 
mann constraint equation (I) becomes 


-2 +y +z? =l. (11) 


The critical points, where x’ = 0, y’ = 0 and z’ = 0, 
correspond to an expanding universe with a scale factor 
a(t) given by a x t?/B(+et)], The effective equation of 
state for the total comic fluid is 


weg = —z? — y’. (12) 


Interaction terms Q have been discussed in the liter- 
ature within the context of inflation and reheating. In 
the conventional reheating model, an interaction term 
T'g¢? dominates at the end of inflation when the scalar 
field is oscillating about the minimum of its potential. 
During this reheating phase the energy transferred from 
the scalar field is completely converted into the matter. 
Within the context of exponential potentials, an interac- 
tion term of the form Q = —ck pm @ was considered in 
Ref. {19}. It was shown that the matter scaling solu- 
tions were stable solutions and the age of the universe is 
older when such the interaction term is included. Certain 
string theories in which the energy sources are separately 
conserved in the Jordan frame naturally lead to interac- 
tion terms in the Einstein frame; scalar tensor_theor 
with matter terms may yield the same results : 
An interaction term of the form Q = 3c H(pp + pm) was 
proposed to look for a dynamical solution to the coinci- 
dence problem fg}. Such the interaction term might be 
motivated by analogy with dissipation. For example, a 
fluid with bulk viscosity may give rise to a term of this 
form in the conservation equation. Without an interac- 
tion term, it was shown that the dark matter could not 
track the phantom energy and would be quickly driven 
to zero id It is of interest to study the cosmological 
consequences of the above two types of interactions in 
the phantom energy model. 


II. MODEL I 


Let us first consider the following interaction 


Q = —CK Pm d, (13) 


where c is a dimensionless parameter. Such a coupling 
arises for instance in string theory ite or after a con- 
formal transformation of Brans-Dicke theory fg}. The 
evolution equations can then be written as an au- 


tonomous system: 
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which has three critical points as follows. 
Point A: 


6 / 2 
va = -Ée ya = 0, ZA = 1+3. (17) 


This solution is physically meaningless since Qm > 1 if 


c#0. 
Point B: 
6 A2 
rp - -Êh yp =j/1+—, zg = 0. (18) 
6 6 
This critical point corresponds to the phantom- 


dominated solution Q, = 1, which always exists for any À 
and c. The effective equation of state, wee = —1 — A? /3, 
depends on the slope of the potential. To find out un- 
der what condition this fixed point is a stable solution, 
we study the behavior of small deviations from the so- 
lution. The linearization of system (4-6) about this 
fixed point yields two eigenvalues mı = — (3 + 7/2) and 


mz = —(àå? +3 — cA). Thus the phantom-dominated 
solution is stable for c < A + 3/.. 
Point C: 
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This fixed point corresponds to the phantom-fluid scaling 
solution, which exists for 


= 2 
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Substituting linear perturbations x —> xc + 6%, y —> 
yc +ôy, and z —> zc +ôz about the critical point into the 
system of equations (4-6), to first-order in the per- 
turbations, gives the following two independent evolution 
equations of the linear perturbations: 
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dy’ = —3yc (5 + zo) ôx — 3y2dy. (21) 


FIG. 1: The convergence of different initial conditions to the 
attractor solution in the (x, y) phase space for the model I 
with c= —2 and à = 1. 


Omega 


FIG. 2: The evolution of the fractional densities of the phan- 
tom field (the dashed curve) and the dark matter (the solid 
curve) for the case in Fig.1. 


The two eigenvalues of the coefficient matrix of the above 
equations determine the stability of the critical point. We 
find that the solution is unstable. 

In the cosmological model with the interaction (I3) 
between phantom field and dark matter, the phantom- 
dominated solution is the only attractor solution in the 
parameter space, c < À + 3/X. In Figs.1-4, we plot the 
numerical results. Comparing Fig.2 to Fig.4, we see that 
the phantom energy more quickly dominates the universe 
when the parameter c decreases. In Fig.1 and Fig.3, the 
trajectories converge at the same fixed point, which is 
only determined by the parameter A. Hence energy trans- 
fer whether from the phantom field to the dark matter or 
vice versa yields the similar cosmological consequences. 
The stable critical point B with £g < 0 indicates that 
the phantom field climbs up the exponent potential. The 
energy density of the phantom field increases as the uni- 
verse expands, which leads to unwanted future singular- 
ity, and therefore the coincidence problem becomes more 
difficult. In the next section, we will investigate a phe- 


FIG. 3: The convergence of different initial conditions to the 
attractor solution in the (x, y) phase space for the model I. 
We choose c= 1 and A= 1. 
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FIG. 4: The evolution of the fractional densities of the phan- 
tom field (the dashed curve) and the dark matter (the solid 
curve) for the case in Fig.3. 


nomenal model, in which the cosmic doomsday is avoided 
and the coincidence problem may be alleviated. 


IV. MODEL II 


Now let us consider the specific interaction 


Q =3cH (pp + pm), (22) 


where c is a dimensionless parameter denoting the trans- 
fer strength. This type of interaction has been proposed 
to look for a dynamical solution to the coincidence prob- 
lem in Ref. bd. Then the equation system can be 


written as an autonomous system: 
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which has four critical points. 
Point A: 
1 
z = 5⁄1 +4c-1), y4 =0, 
1 
2 = 5(v1 + 4c+ 1). (26) 


This solution is physically meaningless since Qp < 0. 
Points B,C,D: The other three critical points are so- 
lutions of the following set of equations: 


f(z) = c, (27) 
y = =g VO, +1, (28) 
z? = 2g? + Eya, (29) 


where we have defined a cubic function 


f(z) =2 (2 + s) ( — Fa.) ; (30) 


The critical point with xg < 0, labelled by B, exists for 


— VA? +12 
24/6 
There are two critical points with xc,p > 0, one of which 


is physically meaningless, labelled by D. The other point, 
labelled by C, exists for 


0<c< f(———— ). 


0<cK< min f(T a"), ne =) 


We see that the point B corresponds to a climbing-up 
phantom field, while the point C corresponds to a rolling- 
down phantom field. In order to study the stability of 
the two critical points, we obtain the two independent 
evolution equations of the linear perturbations 


3 
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FIG. 5: Stability regions of the (A, c) parameter space for 
the model II. In the regions I, either the climbing-up scaling 
solution or the rolling-down scaling solution is the stable late- 
time attractor. In the region II, the climbing-up solution is 
the stable late-time attractor. The solutions are physically 
meaningless in the region III. 


The corresponding eigenvalues of the coefficient matrix 
of the above equations indicate that the critical points B 
and C are always the late-time stable attractor solutions 
if they exist. 

In the case of the interaction form (23), there exist two 
kinds of stable scaling solutions, the climbing-up scaling 
solution with xg < 0 in Figs.6-8 and the rolling-down 
scaling solution with zc > 0 in Figs.9-11. As shown in 
Fig.7 and Fig.10, the universe evolves from the matter- 
dominated phase to the scaling solution, which is charac- 
terized by a constant ratio of the energy densities of the 
dark matter and the phantom field. This may provide 
us with a phenomenological solution of the coincidence 
problem. These results agree with those in Ref. [1 7}. 
The different regions in the (y, c) parameter space lead 
to different qualitative evolutions. In the region II of th 
parameter space in Fig.5, the critical point B is a sta- 
ble solution. However, in the region I both the points 
B and C are stable. Which one is the late-time sta- 
ble attractor solution? The phantom field either climbs 
up or rolls down the exponent potential, which is deter- 
mined by the initial velocity of the phantom field. If the 
phantom field initially climbs up, the effective equation 
of state weg tends to below —1 and realizes a transition 
from Weg > —1 to weg < —1 in Fig.8. Thus the universe 
ends with a big rip. If the phantom field initially rolls 
down, the effective equation of state weg tends to above 
—1 and realizes a transition from weg < —1 to weg > —1 
in Fig.11. In this case the cosmic doomsday is avoided 
and the universe accelerates forever. 


V. CONCLUSIONS AND DISCUSSIONS 


We have presented a phase-space analysis of the evolu- 
tion for a spatially flat FRW universe driven by an inter- 
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FIG. 6: The convergence of different initial conditions to the 
attractor solution in the (x, y) phase space for the model II 
with c = 0.2 and A = 1. We choose initial conditions with 
ro <0. 


FIG. 7: The evolution of the fractional densities of the phan- 
tom field (the dashed curve) and the dark matter (the solid 
curve) for the case in Fig.6. 
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FIG. 8: The evolution of the effective equation of state for 
the case in Fig.6. 


FIG. 9: The convergence of different initial conditions to the 
attractor solution in the (x, y) phase space for the model II 
with c = 0.2 and A = 1. We choose initial conditions with 
zo > 0. 
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FIG. 10: The evolution of the fractional densities of the phan- 
tom field (the dashed curve) and the dark matter (the solid 
curve) for the case in Fig.9. 
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FIG. 11: The evolution of the effective equation of state for 
the case in Fig.9. 
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